A new measure for cosmic shear 



Peter Schneider, Ludovic van Waerbeke, Bhuvnesh Jain & Guido Kruse 

Max-Planck-Institut fiir Astrophysik 
Postfach 1523 
D-85740 Garching, Germany 



Abstract 

Cosmic shear, i.e., the distortion of images of high-redshift galaxies through the tidal 
gravitational field of the large-scale matter distribution in the Universe, offers the op- 
portunity to measure the power spectrum of the cosmic density fiuctuations without 
any reference to the relation of dark matter to luminous tracers. We consider here a 
new statistical measure for cosmic shear, the aperture mass Map(^), which is defined 
as a spatially filtered projected density field and which can be measured directly from 
the image distortions of high-redshift galaxies. By selecting an appropriate spatial filter 
function, the dispersion of the aperture mass is a convolution of the power spectrum 
of the projected density field with a narrow kernel, so that {M^piO)) provides a well 
localized estimate of the power spectrum at wavenumbers s ~ 5/6*. We calculate (Mfp) 
for various cosmological models, using the fully non-linear power spectrum of the cos- 
mic density fiuctuations. The non-linear evolution yields a significant increase of {M^^^ 
relative to the linear growth on scales below ~ 1/2 degree. 

The third-order moment of Map can be used to define a skewness, which is a measure 
of the non-Gaussianity of the density field. We present the first calculation of the skewness 
of the shear in the frame of quasi-linear theory of structure growth. We show that it 
yields a sensitive measure of the cosmological model; in particular, it is independent of 
the normalization of the power spectrum. 

Several practical estimates for (Mfp) are constructed and their dispersions calcu- 
lated. On scales below a few arcminutes, the intrinsic ellipticity distribution of galaxies is 
the dominant source of noise, whereas on larger scales, the cosmic variance becomes the 
most important contribution. We show that measurements of Map in two adjacent aper- 
tures are virtually uncorrelated, which implies that an image with side-length L can yield 
[L/(26')] mutually independent estimates for Map. We show that one square degree of 
a high-quality image is sufficient to detect the cosmic shear with the Map-statistics on 
scales below ~ 10 arcmin, and to estimate its amplitude with an accuracy of ~ 30% on 
scales below ~ 5 arcmin. 



1 Introduction 



Gravitational light deflection caused by an inliomogeneous distribution of matter in the 
Universe causes observable effects on the images of distant sources. Whereas mass con- 
centrations on scales of galaxies and clusters yield strong lensing effects - multiple images, 
(radio) rings, and giant luminous arcs - density inhomogeneities on larger scales or with 
less concentration causes weak lensing effects. In particular, the shape and observable 
flux of distant galaxies is affected by the tidal component of the gravitational fleld and 
the density fluctuations along their lines-of-sight, respectively. Whereas these lensing 
effects are too weak to be detected in individual galaxy images, they can be investigated 
statistically. Assuming that the intrinsic orientations of galaxies are random, cosmic 
shear - i.e., the line-of-sight integrated tidal gravitational field - can be detected as a 
net alignment of galaxy images on a given patch of the sky. 

The statistical properties of the cosmic shear are directly linked to the statistical 
properties of the density inhomogeneities (Gunn 1967, Blandford & Jaroszynski 1981). 
In particular, any two-point statistics, like the two-point correlation function of galaxy 
image ellipticities or the mean quadratic image ellipticity, can be expressed as a redshift 
integral over the power spectrum of the cosmological density fluctuations, weighted by 
geometrical factors depending on the source redshift distribution (Blandford et al. 1991; 
Miralda-Escude 1991; Kaiser 1992, 1996, hereafter K96; Villumsen 1996). These geo- 
metrical factors, as well as the redshift evolution of the power spectrum, depend on the 
cosmological model. Therefore, a quantitative analysis of the cosmic shear statistics can 
provide strong constraints on cosmological parameters and the shape of the fluctuation 
power spectrum. 

It should be stressed that this approach to determine the density fluctuations in the 
Universe does not rely on assumptions about the relation between luminous and dark 
matter. It is therefore of comparable interest as the investigation of the cosmic microwave 
background radiation. In fact, these two approaches are complementary: Whereas the 
CMB measures the fluctuation spectrum at the time of recombination, where the density 
perturbations were small, the cosmic shear probes the fluctuation spectrum at relatively 
small rcdshifts, z ~ 0.5, where it has become non-linear on scales smaller than ~ 5Mpc. 
Secondly, due to the flnite thickness of the recombination shell, the CMB is expected to 
show no structure on angular scales below ~ 5 arcmin, corresponding to comoving scales 
of ~ lOMpc, whereas cosmic shear can probe the density fluctuations also on much 
smaller scales. Finally, a comparison of the results from both techniques can provide a 
beautiful confirmation of the gravitational instability theory of structure growth. 

The above-mentioned authors have calculated the two-point correlation function of 
galaxy image ellipticities and the rms image ellipticity using the linear theory of density 
instabilities. On angular scales below ~ 20 arcmin, one starts to probe the density 
fiuctuation spectrum on scales below ~ S/i"-*^ Mpc (where h is the current value of the 
Hubble constant in units of lOOkm/s/Mpc) where the density fluctuations are nonlinear. 
Jain & Seljak (1997; hereafter JS) generalized the previous investigations of cosmic shear 
to the fully non-linear evolution of the power spectrum, and found that the expected 
rms shear increases by about a factor of two on small scales, relative to the predictions 
from the linear evolution of the power spectrum. Given that the expected rms shear 
on arcminutes scale is of the order of a few percent, this enhancement dramatically 
increases the possibility to detect cosmic shear with currently existing instruments and 
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data analysis techniques. In fact, deep observations of fields around radio QSOs have 
revealed the presence of a coherent shear pattern (Fort et al. 1996, Bower & Smail 1997) 
which has been interpreted as cosmic shear (Schneider et al. 1997), though the selection 
of these fields was based on the magnification bias hypothesis (Bartelmann & Schneider 
1994, Benitez & Martinez-Gonzalez 1997, and references therein) for these luminous radio 
QSOs. Therefore, these measurements should not be interpreted in any statistical sense; 
nonetheless, they have shown that a shear of a few percent is measurable even on scales 
as small as 2 arcmins. 

An attempt to measure cosmic shear on a single ~ 9arcmin field did not yield a 
significant signal (Mould et al. 1994), though a reanalysis of the same data with a some- 
what less conservative approach found a fairly high significance (Villumscn 1995). The 
development of wide-field cameras and sophisticated data analysis methods specifically 
targeted at weak lensing studies (Bonnet & Mellier 1995; Kaiser, Squires & Broadhurst 
1995; Luppino & Kaiser 1996; van Waerbeke et al. 1997) suggest that the discovery of 
cosmic shear in random fields, and its quantitative analysis, are just lurking around the 
corner. 

Higher-order than two-point statistics for cosmic shear are difficult to estimate ana- 
lytically. Bernardeau, van Waerbeke & Mellier (1997; hereafter BvWM) have calculated 
the skewness of the projected surface mass density field, using the quasi-linear theory 
of structure growth. They pointed out that the skewness is a powerful indicator for the 
cosmic density parameter; in particular, it is independent of the normalization of the 
power spectrum, and fairly independent of the cosmological constant. A practical diffi- 
culty related to that measure is that the projected surface mass density is less directly 
linked with observables. Whereas it is tightly related to the magnification of sources, this 
by itself is difficult to observe. One would therefore like to measure the skewness of a 
quantity that is directly related to the shear, which is observable from image ellipticities. 

Such a quantity is provided by the Map-statistics, introduced by Kaiser et al. (1994) 
and Schneider (1996), as a generalization of the C-statistics introduced by Kaiser (1995). 
The latter yields an unbiased estimate of the mean surface mass density within a circle, 
minus the mean surface mass density within an annulus surrounding this circle, and can 
be obtained from the shear within the annulus. First applied to the cluster MS1224 
(Fahlman et al. 1994), this "aperture densitometry" has yielded a lower bound on the 
mass-to-light ratio in the center of this cluster which is considerably larger than values 
typically found in clusters by other means. 

The Map-statistics, which measures the projected density field filtered with a com- 
pensated filter function, combines the properties that it is directly related to the pro- 
jected mass density, and that it is obtainable from the shear, i.e., observable through 
image ellipticities. Furthermore, in contrast to the mean shear within a circle, which is a 
two-component quantity from which no non-trivial third-order moment can be defined, 
the Map-statistics is a scalar whose skewness is well-defined. Schneider (1996) has sug- 
gested that the Map-statistics can be used to search for (dark) matter concentrations on 
high-quality wide-field images. 

In this paper we shall investigate the Map-statistics as a measure for cosmic shear. In 
Sect. 2 we briefiy review the basic equations for the light propagation in an inhomogeneous 
Universe, thereby introducing our notation. The two-point Map-statistics in introduced 
in Sect. 3. It is shown that the dispersion of Map on a certain angular scale can be 
expressed as an integral over the power spectrum of the projected surface mass density. 
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times a filter function containing this scale. This filter function is shown to be quite 
narrow, so that (^M'^^iO)'^ provides a fairly localized estimate of the power spectrum. 
This is contrasted with the rms shear averaged on circles, for which the corresponding 
filter function is very broad. 

For various cosmological models, the dispersion of (Affp) is calculated, for both the 
linear and fully non-linear growth of density perturbations, and compared to the rms 
mean shear. In Sect. 4, we calculate the skewness of Map, following closely the treatment 
of BvWM, i.e., applying the quasi-linear theory of structure growth. Sect. 5 is devoted 
to practical estimators of (M^p) and the skewness, and their respective accuracies. A 
practical advantage of the Map-statistics over the rms shear is that values of Map mea- 
sured in neighboring apertures are virtually independent, whereas the mean shear inside 
a circle has a very large correlation length. This fact is of particular relevance for mea- 
surements of cosmic shear on small angular scales, using wide-field cameras. For a fixed 
total solid angle of available data, the signal-to- noise ratio for the measurement of (M^p) 
is obtained for ^ ~ 2', somewhat dependent on the kurtosis of the density field. Depend- 
ing on cosmology, a measurement of the dispersion of Map with an accuracy of ~ 20% 
should be possible from one square-degree of deep imaging data. We discuss our results 
in Sect. 6. In the Appendix, we consider several approximations which are used through- 
out the main body of the paper. We focus on the contributions to the skewness which 
arise even in the case of Gaussian density fluctuations, and find that by including these 
effects, the skewness changes by ~ 10%, which justifies the use of our approximations. 



2 Light propagation in slightly inhomogeneous Universes 

We shall use a notation similar to K96 and JS. The metric of the homogeneous Universe 
is written in the form 

ds^ = dt" - a''{t) [dw^ + fl{w)du;^] , (2.1) 

where a{t) = (1-1-2;)""'^ is the dimensionless cosmic scale factor, normalized to unity 
today; is a radial coordinate, and fxiw) is the comoving angular diameter distance to 
distance w. The spatial curvature 

K = (f2d + , (2.2) 

is related to the present day value of the density parameters and in dust and in 
vacuum energy, in terms of which fxiw) reads 

sin (yKu?j for > , 

w for K = , (2-3) 

(-K)-i/2 sinh (a/^w) for i^: < . 

The angular diameter distance D{zi, Z2) is conveniently expressed in terms of fx by 

D{zu Z2) = j^—fKiw2 -wi) , (2.4) 
zi < Z2, where the Wi are the radial distances corresponding to redshifts Zi, 
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^= rdz'^=^r ^ (2.5) 

Jo H Ho Jo V(l + ^')^^d + (1 + -2')^ (1 - - ^v) + 

Light propagation in a weakly inhomogeneous Universe has been investigated in many pa- 
pers (e.g., Blandford et al. 1991; Seitz, Schneider & Ehlcrs 1994, and references therein). 
We shall use the results of these papers in the following form: 

Consider a bundle of light rays intersecting at the observer. Each of these rays 
is characterized by the angle d it encloses with a fiducial ray. Let x(0, w) denote the 
comoving transverse separation of the ray characterized by from the fiducial ray; it 
satisfies the propagation equation 



(x (6>, k;) , w) - V±^(°H«;)1 , (2.6) 



where (x (0, vS) , w) is the Newtonian gravitational potential at comoving distance w 
and comoving perpendicular separation x (0, w) from the fiducial ray, <?(°) {w) is the po- 
tential along the fiducial ray, and V_l = ^gf^, '^l) transverse gradient operator 
in comoving coordinates. Here we have assumed that the fiducial ray propagates nearly 
parallel to the local ,T3-direction, and since all angles involved are small, all rays con- 
sidered propagate nearly parallel to this direction. If the Newtonian potential vanishes, 
x(0, ly) ~ fK{w)d, in agreement with the identification of fxiw) as comoving angular 
diameter distance. A formal solution to (2.6) is obtained from the Greens function of 
the operator on the left-hand side of (2.6), yielding 

2 /•"-' r -1 

x{0,w) = fK{w)e-^ dw' fK{w-w') V±${x{0,w'),w')-V±^^^\w') .(2.7) 

c Jo L J 

A source at w with comoving distance x from the fiducial ray will be seen in the absence 
of light defiection at an angle (3 = x./ fK{w), which we shall call the unlensed source 
position. Defining, as in usual lens theory, the Jacobian matrix A as 

..^ X 93 1 9x , , 

^^'•^^ = W = 7^)80 • <^-«> 

we obtain from differentiation of (2.7), 

Aij {0, w) = r dw' fK(w-w')fKiw') ^^^^ ^ ^ ^2.9) 

c Jo Jk[w) 

where indices z on ^ preceded by a comma denote partial derivatives w.r.t. Xi. In general 
this is not an explicit equation for ^, since in order to calculate A, one first has to solve 
for the ray position x(0, w) and then solve the integral equation (2.9) for A. However, 
for weak gravitational fields which are of interest to us here, one can expand A in powers 
of the Newtonian potential and keep only the lowest order term; this results in 

A,{0M = - \ r^^' ^ (/^(«;')0,^') . (2.10) 

Jo fK{w) 

Hence, to linear order in the distortion is obtained by integrating the second derivatives 
of the potential along the unperturbed ray. This feature allows for a dramatic simplifi- 
cation of the calculations below. In the Appendix, we shall consider some higher-order 
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terms for A corresponding to lens-lens coupling and to dropping the 'Born approxima- 
tion'; i.e., we shall calculate A up to second order in Also note that A given in (2.10) 
is symmetric. One can therefore define a defiection potential, 

m^) = |/"d»' M^i^^fuM^') . (2.11) 

in terms of which one can treat lensing by large-scale structures similarly to the single 
lens-plane case; e.g., the corresponding 'surface mass density' k{0,w) and shear 7(0) = 
71 + 172 are given as 

k{9,w) = - (V'.ii + V',22) , 

J (2.12) 

where indices on ifj separated by a comma denote partial derivatives with respect to 9i. 
In particular, Aij = 6ij — Thus, one can consider the surface mass density k - or 
its associated defiection potential ip - as the fundamental quantity, and in fact it is the 
only quantity which is probed by cosmic shear measurements. One obtains from (2.11) 
and (2.12), by adding a term 33 to the integrand which cancels out upon lu-integration, 
and by using Poisson's equation in the form 

V^^ = ^I^Mks , (2.13) 
the following expression for the projected density field: 

2 V c y Jo fK[w) a{w') 

where 5(x, w) is the density contrast. The projected density field depends on the source 
redshift (or distance). We shall assume that one observes the shear through a population 
of galaxies for which only the redshift probability distribution pz (z) is known, or equiv- 
alently, pyj{w)dw = pz{z)dz. Then, the source distance-averaged projected mass density 
becomes 



dw pw{w) k{0,w) = - i-^j Od J dwg{w)fK{w) 



SifK{w)d,w) 



a{w) 

(2.15) 
where 

g{w) := r dw' pUw') ^^^^^^ (2.16) 
Jw Jk{w') 

is the source-averaged distance ratio Dds/Dg for a density fluctuation at distance w, and 
wu is the comoving distance to the horizon, obtained by putting 2; = cxd in (2.5). In this 
paper, we shall consider two different redshift distributions of sources. In the first case, 
all sources are assumed to reside at the same redshift Zg, so that 

Pz{z) = 5T)iz - zs) . (2.17) 

More realistically, we consider a redshift distribution of the form (e.g.. Small et al. 1995) 
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Pz{z) oc z'^ exp {-[z/zof) 



(2.18) 



The mean redshift of this distribution is proportional to zq and depends on the parameter 
P which describes how quickly the distribution falls off towards higher redshifts. In 
particular, for /5 = 1.5 (a value that we shall use throughout), (z) = I.SOSzq. 

3 Two-point Map-statistics 

3.1 The power spectrum of the projected density field 

Provided the density contrast 5 is a homogeneous and isotropic random field, so is the 
projected density k. Consider the Fourier transform of the projected density field, 

R{s) := j d'^e K{e)e-^^-^ . (3.1) 

We define the power spectrum Pk,{s) of k by 

(«(s)«*(sO) = (27r)2<5D(s - s') P,(|s|) , (3.2) 

where the Dirac delta 'function' 5d occurs because K,{d) is a homogeneous random field, 
and the dependence of on the modulus of s only expresses the fact that k, is an 
isotropic random field. Following K96, we can calculate Pk{s) from the power spectrum 
of the density fiuctuations, defined accordingly by 

(S{k)6*{k')^ = {2nf6B{k-k')P{\k\) . (3.3) 

The derivation of the Fourier-space analog of Limber's equation in K96 is valid provided 
the power spectrum P does not evolve appreciably on timescales corresponding to the 
light travel time across the largest significant fiuctuations, and provided the typical source 
distance is much larger than the largest-scale fiuctuations. With these two assumptions, 
K96 (see also Kaiser 1992) obtains 

P.(.) = ?(M^«.r«d„4Mp(^,..). (3.4) 

4 V c y Jo a^{w) \fK{w) J 

The second argument of P indicates that the power spectrum evolves (slowly) with 
redshift. We shall later derive an analogous relation for the three-point function, using 
the same strategy as in the derivation of (3.4). Since s is the Fourier-conjugate of the 
angle 9, we can relate an angular scale to its corresponding s by s = 27r/^ = 2.16 x 
10^ (6'/larcmin)-^ 

In Fig. 1, we have plotted Pk{s) for five different cosmological models. For three of 
them, the power spectrum P{k) is approximately cluster normalized, which corresponds 
to a"8 ~ 0.6 for an Einstein-de Sitter universe (EdS, — 1, — 0), a% — 1 for 
both an open universe (OCDM, = 0.3, i?v = 0) and a spatially fiat universe with 
cosmological constant (ylCDM, f2d = 0.3, f}^ = 0.7). In all these cases, we have used 
the CDM spectrum as given by Bardeen et al. (1986), but set the shape parameter of 
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the linear power spectrum to F — 0.25, which yields the best fit to the observed two- 
point correlation function of galaxies (Efstathiou 1996). The remaining two cosmological 
models have EdS geometry, but either a higher normalization (erg = 1, approximately 
corresponding to the COBE normalization) or a different shape parameter — 0.5, 
corresponding to the original definition for T if i^o = 50km/s/Mpc). For each model, 
the projected power spectrum has been calculated for a linearly-evolved cosmological 
power-spectrum (thin curves) , and for the fully non-linear power spectrum, following the 
prescriptions of Hamilton et al. (1991), Jain, Mo & White (1995), and Peacock & Dodds 
(1996) - we used the fit formulae of the latter paper throughout. 




10-2 0.1 1 101 IQ2 103 10" 10' 102 103 104 105 106 
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Fig. 1. The power spectrum Pk{s) of the projected density field, as defined in (3.4), for five different 
cosmological models, as indicated by the line types; the numbers in parenthesis are (as,r). For each 
cosmological model, the thin curves correspond to the linear evolution of the power spectrum P{k,a), 
whereas the thick curves are calculated with the fully non-linear evolution of the power spectrum, as 
given in Peacock & Dodds (1996). For this plot, the redshift distribution (2.18) was used, with zq — 1 
and f3 = 1.5. Note that the large amplitude of the EdS model with 7 = 0.25 relative to the OCDM and 
ylCDM model is due to the factor in (3.4) 



3.2 The Map-statistics 

Define the aperture mass by 

Map(^):= J d^i9U{\^\)K{^), (3.5) 

where the integral extends over a circle of angular radius 6, and U{'d) is a continuous 
weight function which vanishes for d > 6. Provided t/ is a compensated filter function, 
i.e., 
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d'd'dU{'d)^0 , (3.6) 







one can express Map in terms of the tangential component 7t of the shear inside the 
circle (Kaiser et al. 1994; Schneider 1996), 

Map(^)= J d2^Q(|^|)7t(t?), (3.7) 

where 

2 /"^ 

<5(^) = ^y d'd' 'd'U{^')-Ui^) , (3.8) 

and the tangential component of the shear at a position = (?? cos (p, sin (p) is 

7t(^) = -7^e(7(^)e-2i^) . (3.9) 

The useful property of Map is thus that, one the one hand, it yields a spatially filtered 
version of the projected density field, and on the other hand, that it can be expressed 
simply in terms of the shear. Since in the weak lensing regime, the observed galaxy 
ellipticities provide an unbiased estimate of the local shear. Map is directly related to 
observables. 

Obviously, the ensemble average of Map vanishes, (Map) = 0. The dispersion of Map 
can be calculated as follows: 

(Mfp(^)> = / u{e') j dH u{^) («(0')«W> 

= / dH' UiO') J d^^ Ui^) J ^ e-(^'-^)p.(.) (3_io) 



= 27r ^ ds s P^{s) d-d ^ Ui-d) Jo(s'i?) j , 



where we used in the first step that the two-point correlation function is the Fourier 
transform of the power spectrum, and the Bessel function Jq comes from the angular 
integrations in 0' and t?. 

The choice of the weight function U{'d) is arbitrary at this point. We shall write 
U{'d) = u{'d/9)/e^, with u{x) = for a; > 1, and Q(??) = q{'d/9)/e^. Furthermore, we 
choose the normalization such that 

271 / dddQ{d)^2TT ! dxxq{x) = l. (3.11) 
Jo Jo 

A set of weight functions which satisfy (3.6) and (3.11) is 

(£ + 2) 



u{x) 



which peak at a; = and go to zero with order £ as a; — > 1. Correspondingly, 

TT 
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Then defining 



= ^ -V ' h+M , 

TT 

we can write the dispersion as 

/•OO 

(Mfp(^)> = 27r / ds s P^is) [leisd)]' 
Jo 



(3.14) 



(3.15) 
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Fig. 2. The filter function lj{r]) in the dispersion of Map for three different value of the t. I = \ (dashed 
curve), 1 — 2 (dotted curve), and = 3 (solid curve). For comparison, the corresponding filter function 
for top-hat filtering /.|jj(?7) is also plotted (dashed-dotted curve). For larger values of 77, all functions 
rapidly oscillate, owing to the Bessel functions. For clarity, we have therefore plotted the amplitude of 
these oscillations as corresponding thick curves 



We have plotted the filter function /|(r7) in Fig. 2, for three different values of i. 
One can see from (3.14) that Ii{r\) oc for small r/, and that for large r/, the function 
Ii{r\) oscillates with an amplitude oc Hence, the filter function /| behaves like 

'xf and in the two respective hmits, and is indeed a very locahzed filter, with a 

peak at ?7 ~ 5 for small t. Hence, {M1^{Q^^ provides an accurate measure of the power 
spectrum of the projected density field at s ~ 5/6*, with little dependence on the local 
shape of this power spectrum. 

We shall compare (Mfp(6')) with the rms value of the shear averaged over a circular 
region of angular radius d\ let 
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7(^):=^ J d^^7W, (3.16) 

then its dispersion is 

J7|' (^)) = 271 dssP,{s)[lTii{se)]\ (3.17) 
with 

/TH(r7) = ^ (3.18) 

777] 

(cf. Blandford et al. 1991). The dash-dotted curve in Fig. 2 displays for comparison 
with l£. One sees that /^jj is a much broader function which tends towards a constant for 
rj 0, and its amphtude decreases hke r]~^ for large ij. Therefore, the shear dispersion 

7|^^ is a much coarser probe of the power spectrum than (M^p(6')). 
3.3 Results 

In Fig. 3 we have plotted the rms values of the Map-statistics (upper panel), as well as 
that of the mean shear within circles (lower panel). The same cosmological models as 
in Fig. 1 were used, and we present results both for the linearly evolved cosmic power 
spectrum (thin curves) and for the fully non-linear evolution (thick curves). Whereas 

/ 2\^^'^ 1/2 

(|7| ) decreases monotonically with 0, the shape of (M^p) closely reflects the shape 

of the projected power spectrum displayed in Fig. 1; this is of course related to the 
narrowness of the filter Ii shown in Fig. 2. In Fig. 3, and for the figures of the remainder 
of the paper, we used £ = I, but have checked that changing to £ = 2 or 3 does not yield 
qualitatively different results. Compared to the prediction of the linear evolution of the 

1 /2 

power spectrum, the peak of (M^p) is shifted to substantially smaller angles; at the 

1 /2 / 2\^^'^ 

same time, the non- linear evolution affects (M^p) more than ( |7| ) , as the latter 
picks up power from larger scales which are hardly affected by non-linear evolution. 

1/2 

As can be seen from the figure, the values of (Affp) are substantially smaller than 

/ - 2\^/2 

those of ( I7I ) , at least on scales below one degree. This is related to the fact that at 



a given angular scale 9, the Map-statistics is sensitive to smaller-scale structures than the 
|7|-statistics, as can be seen in Fig. 2, and so these two rms values should be compared 
at different 'effective' scales. Whereas this difference in magnitude may suggest that 
the Map-statistics is observationally disfavoured, we shall show in Sect. 5 that it has the 
advantage that measurements in neighboring fields are nearly uncorrelated. 

The dependence of the rms of Map on the source redshift distribution is displayed 
in Fig. 4, for the same five cosmological models, and two angular scales. The dependence 
on Zs is generally weak, in rough agreement with the power law dependence z^ -^ found by 
JS for the top-hat filter. The model with cosmological constant shows a stronger growth 
with ^s, as distances grow faster with redshift in such models, and thus provide larger 
path lengths for lensing. 



4 Three-point A^ap-statistics: skewness 
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Fig. 3. The rms values of the shear are shown vs. filter scale for the same cosmological models as used 
in Fig. 1. The upper panel shows the rms value of Map with the filter Ii defined in (3.14), while the 
lower panel shows the rms shear computed with a top-hat filter (3.18). Note the different scales in the 
two panels. The source galaxies are assumed to follow the distribution (2.18), with (3 — 1.5 and zq — 1. 
The thin curves display the prediction for the rms shear if linear evolution of the density fluctuation 
spectrum is assumed, whereas the thick curves follow from the fully non-linear evolution of the power 
spectrum. The results shown in the lower panel are fully equivalent to those of JS, except that they 
considered a redshift distribution of the form (2.17), and they considered the 'polarization' as a measure 
of net galaxy ellipticity, which in the weak lensing case equals twice the shear 



If K is a Gaussian random field, so is Map, and in particular, the expectation value of the 
third-order moment of Map would vanish. However, since the nonlinear gravitational evo- 
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Fig. 4. The rms value of Map as a function of source redshift, for the same five cosmological models as 
in Fig. 1, indicated in the figure; the numbers in parenthesis are (erg, F). The thick curves are for 9 = 1', 
the thin curves for = 10'. The fully non-linear evolution of the power spectrum has been used. The 
sources were assumed to be all at the same redshift 

lution of density fluctuations transforms an initially Gaussian field into a non-Gaussian 
one, this third-order moment is non-zero in general. In fact, (Mfp) is a measure of the 
non-Gaussianity of the density fluctuations at medium redshifts. Since 



(Mfp(^)> = J d^^i u{ei) J d^e^ u{e2) J d^e^ u{e^) (ac(0i)av(02)«(03)) 

= j d^Oi u{ei) J d2^2 u{e2) J d^Os u{9^) (4.1) 



(27r)2 J (27r)2 J {2n) 

the evaluation of the third moment requires the calculation of the three-point correlation 
function of k. 

The skewness of the projected surface mass density has already been discussed in 
BvWM. They have calculated the skewness of the projected density field by considering 
the quasi-linear evolution of the density fluctuations. Below, we shall follow the same 
procedure to calculate (Mfp). As was also pointed out by BvWM, even for the strictly 
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linear evolution of the power spectrum, when the density field conserves its initial Gaus- 
sian nature, the observable skewness would not vanish identically. We get back to this 
point at the end of this section. 



4.1 The three-point correlator for k 

We shall now evaluate this three-point function in terms of the corresponding function 
of the cosmic density field, i.e., to obtain the analog of (3.4) for the three-point function. 
For abbreviation, we write (2.14) in the form 

K{e,w)= dw' G{w,w')5{fK{w')d,w') , (4.2) 
Jo 

with 

^ (Ho\\ fK{w-w')fK{w') 

G{w,w) = - — i2a- 



^\cj fK{w)a{w') 
for w' < w, and zero otherwise. Then, 

(k(si, u;i)k(s2, W2)k{s3, Ws)) 

= Jd^O^ e-^^-^^ Jd^02 e-i^-^^ Jd% e-^^-^3 



X 



/•wi /•W2 rw-s 

/ dVi G{wi,Vi) / dV2 G{W2,V2) / dVs G{W3,V3) 

Jo Jo Jo 



(4.3) 



X {S{fK{vi)0uVi)5{fK{v2)02,V2)S(fK{v3)es,Vs)) , 



where we have inserted the Fourier transforms of ii and used (4.2). Assuming that the 
largest scale fluctuations are much smaller than the typical distance to a source, the 
three-point correlation function of S will vanish except when vi ^ V2 ^ v^. Thus over 

the t^-range where (SSS) does not vanish, we can set the t^-arguments in G to be all 
the same, and we can also set all fK{vi) equal. Since G{wjv) = if v > w, the outer 
^-integration only extends to Wmin = m.m{wi,W2,ws). Hence, after replacing d by its 
Fourier transform, one obtains 

{k{si, wi)k{s2, W2)k{ss, ws)) 

= J d^e, e-'s-^^ J d^e2 e-'s-^^ J d% e-'^-^3 

dv G{WI,V)G{W2,V)G{W3,V) J dv' J dv" 

r d^kl i/^(,;)0i.ki ifcist; f d^^2 ifK{v)02-'k2jk23v' 

J (27r)3 J (27r)3 

/ §^ e^^-(-)^-k33ife33." (^S{k^)~^{k2)~^{k3)) , 

where we have written the 3-dimensional vector k as (k, ks). The 6'i-integrations can 
now be carried out, each yielding {27r)'^5r){si — fKiv)\s.i). After that, the kj-integrations 
become trivial. The v' and t'"-integrations can be carried out yielding delta-functions in 
k23 and k^s, which are trivially integrated away. One thus finds: 



(4.4) 
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{k{Si, Wi)k{s2, W2)k{ss, W3)) 



dv G{wi, v)G{w2, v)G{w3, v) 



X 



1K3V 



(27r) 



k. 5 



S2 



s 



S3 







(4.5) 



Next, the average of (4.5) over a source redshift distribution is performed. The only point 
to notice here is the upper hmit of integration for v. Integrating (4.5) over Y[ J Pw{wi) 
(recaU that Pwiw) dw — Pz{z) dz is the source redshift distribution) gives an expression 
of the form 



/■Wh rUin rWn r 

= / dwi / dw2 / dw3 

Jo Jo Jo Jo 



dv F 



dv 



dwi 



dW2 







dw3 F . 



In obtaining the second equality we have made use of the fact that the integrand F is 
symmetric in wi,W2,ws, and that G{w,v) = if v > w. We therefore obtain for the 
redshift-averaged three-point correlation 



(k(Si)k(S2)k(S3)) = Y 



X 



dks 



(27r) 

with g{w) as defined in (2.16) 



Sl 



/k(w) 





,k3 



dw 
S 



a^{w)fl^{w) 

S2 



fKiw) 



,0 s 



S3 



fKiw) 



,0 



(4.6) 



4.2 Quasi-linear theory of density fluctuations 

In order to calculate the triple correlator in (4.6), we shall use quasi-linear perturbation 
theory, in which the density field 6 is considered as a 'small' quantity and expanded into 
a perturbation series, 6 = 6^^^ + 5*^^-* + . . ., where 5'-"-' = O ([5*-^'*] )• Here, (5*-^-* is the 

linearly evolved density perturbation, d^^\k,w) = D^{w)SQ^\k), D^{w) is the linear 

growth factor, normalized to -D+(0) = 1, and SQ^\k) is the density perturbation linearly 
extrapolated to the present epoch. The perturbation series can then be inserted into the 
continuity equation and the Euler equation, and a closed solution for every order can 
be obtained in terms of lower-order terms (see, e.g.. Fry 1984, Goroff et al. 1986, and 
references therein). In particular, in an EdS Universe, for the first non- linear term, one 
obtains 



~S^'\k,w) = Dliw) I ^ ~sl,'\k') ~dl,'\k-k') 



X 



5 

7 + 



(27r)3 
k' -ik-k' 



k' 



2 

+ 7 



k' ■ k-k' 
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k' 



k-k' 



(4.7) 



— * 

Bouchet et al. (1992) showed that the ^-dependence of this term depends weakly on 
cosmology, and that it is an excellent approximation to restrict the cosmology dependence 

solely to the growth factor D\{w). In order to calculate (sss"^, we note that the tri- 

linear correlation of the linear density field vanishes, owing to its assumed Gaussian 
nature. Thus to lowest order, we have 
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+ 2 terms obtained from permutation 



(4.8) 



Considering only the first term, inserting (4.7), making use of the fact that 5^^\k,w) 
D^{w)SQ'\k), and using the relation 



+ 2 terms obtained from permutation 
= {27Tf (^Po{ki)Po{k3)5u{h + k2)Su{k3 + ^4)) + 2 terms , 
valid for Gaussian fields, we find 

~5^^\h)~5^^\k2)5^^\h)) (w) = 2{27rfDX{w)Po{k,)Po{k2)5-D{h + ^2 + ^3) 



(4.9) 



X 



5 1 
7 + 2 



v 



ki 



ki ■k2 + 



ki ■ k2 



1 2 n 



ki 



(4.10) 



Here, Po{k) is the power spectrum of the Hnearly extrapolated density field. This expres- 
sion is only one of three terms appearing in (4.8), with the other two being obtained by 
permutation of the k^. In order to calculate M^^ according to (4.1), we have to integrate 
over all s^, and each of the three terms in (4.8) yields the same contribution. Therefore, 
we can just use (4.10) instead of (4.8), and multiply the result by a factor 3. Then, 
combining (4.1), (4.6) and (4.10), we obtain 



(M.„W> = -(2.) ^-j «,y^ A,n „3(„);^(„) 



X /(|Si +82!^) 



5 1 

7 + 2 



1 1 



Si 



S2 



, 2 [si • S2] 
Si ■ S2 + -- — r^r. — 



Pl S2 



4.3 Skewness 

We define the skewness as 

We have calculated S{9) using the results of the preceding subsection. Quasi-linear theory 
is expected to underestimate (Mfp(^)), in particular on scales below ~ 20 arcmin which 
roughly demarcates the linear and nonlinear regimes. However it has been shown that the 
quasi-linear theory yields a surprisingly accurate estimate of the skewness of the density 
field, even in the regime {\S\ ~ 1) where one could not expect the perturbation series to 
yield reliable results (e.g. Baugh, Gaztaiiaga & Efstathiou 1995, Colombi, Bouchet & 
Hernquist 1996, Gaztafiaga & Bernardeau 1997). Therefore, we expect that the results 
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for S{9) as calculated here will provide a good approximation to the true skewness, even 
on scales down to a few arcminutes. On very small scales with density contrasts much 
larger than 1, the skewness of the density measured in N-body simulations is indeed 
larger than the perturbation theory value, but by no more than a factor of two (Baugh 
et al. 1995, Colombi et al. 1996). 
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Fig. 5. The skewness, as defined in (4.12), as a function of filter scale 6, for four cosmological models. 
Since the skewness as calculated from quasi-linear theory is independent of the normalization of the 
power spectrum, the two EdS models with the same 7^ yield the same S(0) curves. The sources were 
distributed according to (2.18), with 2:0 = 1 and /3 — 1.5, and the filter with £ — 1 was used 



In Fig. 5 we have plotted the skewness as a function of angular scale, for various 
cosmological models. Going from an EdS model to an open model, the skewness increases 
by more than a factor of two, with the spatially flat A model taking intermediate values. 
The two EdS models with different power spectra have only slightly different skewness. 
This result is in full agreement with BvWM, who found that the skewness is a fairly 
sensitive probe for the cosmological model, to a large degree independent of the exact 
shape of the power spectrum. The skewness on an angular scale of 5' is plotted in Fig. 6, 
for cosmological models with i?^ = 0, and flat cosmologies (i7v = 1 — ■^^d); and two 
different shape parameters F of the power spectrum, as a function of i7d- The two panels 
are for different source redshift distributions. Again we see that the variation of S with 
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the cosmological parameters is much stronger than with the shape of the power spectrum. 
And by definition S is independent of the normahzation erg of the power spectrum, at 

least in perturbation theory, where it is the ratio of two terms, each of O (^[^^^'^]'^^ ■ 

Wc therefore conclude that the skewness is a powerful discriminator between different 
cosmological models. 

Since (Affp) contains a factor i7|, and (Affp) has a factor one expects that to 
leading order, S oc i^^^- As can be seen from Fig. 6, S follows this expectation rather 
closely for the models with = 0. For flat cosmologies, however, these f2d factors are 
no longer the dominant terms in the dependence of S on f2d- The dependence on is 
weakened due to the distance factors in the integrands for (Affp) and (Mfp). 

Up to now, we have considered the skewness as it arises due to the quasi-linear 
evolution of the density field which transforms an initially Gaussian field into a non- 
Gaussian one. The linear density field would not cause any skewness in the frame of 
the approximations used up to now. Nevertheless, even in the case of Gaussian density 
fiuctuations, the observable skewness would not be identically zero, owing to the following 
three effects: (i) Light rays do not propagate along 'straight' lines, but are deflected. 
Therefore, the separation vector x{6,w) deviates from fK{w)0. This deviation from 
the 'Born approximation' leads to a non-vanishing skewness. (ii) In the transition from 
(2.9) to (2.10), in addition to the Born approximation the matrix A on the right-hand- 
side of (2.9) was approximated by the unit matrix; in that way, the leading order term 
(in of the propagation matrix was obtained. The next order term, which describes 
the coupling between lens planes, yields a non-vanishing skewness. These two effects 
have been quantitatively analyzed in BvWB for the skewness of the projected density 
fleld, where it was shown that they are small compared to the expression obtained from 
quasi-linear density evolution, (iii) The fact that the image ellipticity is not an unbiased 
estimator for the shear 7, but for the reduced shear g = 7/(1 — «;), implies that a Gaussian 
density fleld can produce a flnite skewness to the extent that the weak lensing condition 
K <C 1 is no longer valid, and that we have to consider the term We shall 

consider all three effects in the Appendix. All three effects yield contributions to (Mfp) 
which are proportional to the square of the power spectrum Po{k), as is the case for 
the leading term considered here. Therefore, there is no a priori reason to expect that 
these 'correction' terms are much smaller than (4.11). As it turns out, however, these 
correction terms amount to ~ 5% of the leading term (4.11), and may therefore be safely 
neglected presently. 

A flnal point concerns the redshift distribution of the sources. Fig. 6 shows that 
the source redshift dependence of S is rather large, as pointed out by BvWM. This is a 
potential problem for two reasons: (1) First, it means that in order to compare an ob- 
served skewness with theoretical calculations, we have to know the redshift distribution 
fairly precisely. Since the likely sources are very faint, a complete spectroscopic survey 
to the corresponding magnitude limits is not available, but progress in the method of 
photometric redshifts may yield sufficiently accurate estimates of the required distribu- 
tion. (2) The magniflcation bias changes the redshift distribution of the observed galaxy 
iamges in a way which depends on the local value of the magniflcation. The observed 
skewness is sensitive to this. It is not the aim of this paper to discuss these points. A 
rough estimate of the second point shows that the magniflcation bias is negligible, but a 
detailed and systematic analysis of these effects remains to be done. 
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Fig. 6. The skewness, as defined in (4.12), as a function of i?^, for vanishing cosmological constant 
i7v = 0, and for fiat Universes, i7v = 1 — ^d- For both cases, two values of the shape parameter F were 
considered. The redshift distribution of the sources was assumed to follow (2.18), with f3 = 1.5, and 
ZQ — 0.65 (upper panel) and zq = 1.3 (lower panel), corresponding to a mean redshift of about one and 
two, respectively. As can be seen, the dependence of S on the cosmological parameters is considerably 
stronger than the dependence on the shape of the power spectrum 

5 Practical estimates 

We shall now consider some practical estimators for the Map-statistics, and their dis- 
persion. The results of this section yield an estimate of the observations necessary to 
measure (M^p) with given precision. 
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5.1 Dispersion for a single field 

Consider first a single circular field of angular radius 9, in which N galaxies are observed 
at positions Oi with complex ellipticity e^. In the case of weak lensing, k <^ 1, the 
transformation from the intrinsic ellipticity e^^^ to the observed one is simply e = 7 + e*^^-* . 
Since the intrinsic orientation of galaxies is random, e is an unbiased estimator of the 
local shear. Defining in analogy to (3.9) the tangential component of the ellipticity of a 
galaxy at di = {9i cos (fi, di sin (^j) by 

eti = -7ee(eie-2i^') , (5.1) 

then eti is an unbiased estimator of 7t(0i) = 7ti- In terms of the observable image 
ellipticities, we define the estimator 

where Qi = Q{di), which in turn is defined in (3.8), and the sum is taken only over 
terms with i ^ j. The expectation value of M is obtained by three averaging processes: 
averaging over the intrinsic ellipticity distribution, averaging over the galaxy positions, 
and the ensemble average. We denote the first of these processes by the operator A, and 
the second by P; these two operators commute. Then, 

2 

A(et,et,-) = A ( [7^, + e^^^] [7ti + eg^] ) = lultj + ^^ii , (5.3) 

since A(e^^'*) = 0, owing to the random intrinsic orientation; here, (jg is the dispersion of 
the intrinsic ellipticity distribution. The term with is divided by 2 as we want the 
dispersion of only one component of the complex ellipticity. In (5.3) we also used the 
fact that the intrinsic ellipticity is uncorrelated with the shear. Note that the second 
term in (5.3) does not appear in the sum of (5.2). The operator P which averages over 
galaxy positions is defined as 




P(^) =\\ \ i (5.4) 
so that 

P I QiQ^lultA = ^^7^^ J d'^i Qi^i) J Q(^2)7t(0i)7t(02) , (5.5) 

because the sum yields A'"(A'" — 1) equal terms. The expectation value E(M) of M then 
becomes: 

E(M) = (P(A(M))) = J d'e, Q{e,) J d'd^ Qie^) {^M^M) = . (5.6) 

Hence, M is an unbiased estimator for (Mfp). 

Next, we consider the dispersion of this estimator. 
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(7\M) = E(M2) - [E(M)] , (5.7) 

with 

= j^^L _ E Q^Qj^tietj Yl QkQietketi • (5.8) 
Starting with performing the average over the intrinsic eUipticity distribution, we find 

./ \ , « (s) (s) (s) (s)^ 

A(etietjetfcetO = ^ultjltkltl + A [el/ el/ el,^ el/ j 

(5.9) 

+ Y i^^i^^j^kl + lUltk^jl + Itiltl^jk + ItjltkSil + ItjltAk + ItkltAj) ■ 

Owing to the restrictions in the sums of (5.8), the first and last term in the second 
parenthesis do not contribute. 

To evaluate the position-average over the first term in (5.9), we have to consider 
three difi'erent cases: (a) all four indices i, j, k, I are different; (b) one of the second pair 
of indices k, I is equal to one of the first pair (c) both indices in the second pair are 
equal to those in the first. These cases occur in the sum of (5.8) in A^(A^'— 1)(A^'— 2)(A^'— 3), 
4:N{N - 1){N - 2), and 2N{N - 1) terms, respectively. Therefore, 



N N 



] 



p I 5Z X] QiQjQkQiitiitjiikiti 

yi,3i^i k,l^k j 

A-(A--l)(A--2)(A--3) 4 4A-(A--l)(A--2) 

(7r^2)4 ap + ^^^2)3 

2 



Mfp/ d¥g2(^)7t'(t?) 

(5.10) 

For the second term in (5.9), we note that since j^k ^l, k (^^u e^/^ e/^^ e/^^ contributes 

to the sum in (5.8) only for k = i and I = j, or k = j and I = i. There are 2N{N — 1) 
such terms. Hence, 

p (t t Q'Q.Q^Q''^ (4f4?4*'4?)) = (f )'«^ ■ (8-") 

y^ii?^* k,l^k J 

where we have defined 

G = 7Te^ J d^^Q^i^) . (5.12) 

Finally, each of the four terms in the last term of (5.9) which contribute to the sum in 
(5.8) yield the same result when averaged over position. Considering one of these, 

(N N \ f ^ 

Y Y QiQjQkQiitiitk^ji = P X] QiQ]Qkitiitk 
i,j^ik,ij^k J \j,i^j,k^j ) (5.13) 

_iV(Ar-l)(Ar-2) ^ . ^ , Ar(Ar-l) 



(7r^2)< 
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where the first term comes from terms with i ^ k, and the second are those from i = k. 
CoUecting terms, and taking the ensemble average yields for the dispersion 

^ (N-2)iN-Z) m-2) 

^ ' N{N - 1) \ N{N - 1) ^ ' / 

2 o / - 1 



with 



N{N - 1) ' \ « / ^ 2N{N - 1) ' \ ap/ > 

Note that the dispersion contains three different sources of noise: the contribution due to 
the finite width of the intrinsic ellipticity distribution, the 'cosmic variance', and the noise 
due to the finite number of randomly located galaxy images. The latter effect is contained 
in the terms which include Mg. However, this source of noise is never dominant: since 
the magnitude of will be comparable with that of Mfp, the second term in (5.14) 
is smaller by a factor 1/N than the first, and smaller by a factor 2GM^/a1 than the 
fourth term. Similar estimates are valid for the other two terms containing M^. Thus, 
by dropping terms containing Mg and taking A?" 1, we find 

a^(M)«.^4(Mi>V(^ + x/2(Mfp>y , (5.16) 



where /X4 is the kurtosis of Map, //4 = (-^ap) / (-^ap) ~ ^1 which vanishes if Map is 
distributed like a Gaussian. We note that G is a factor of order unity; if we choose the 
weight function (3.12 & 13), we obtain 

(l + £)(2 + £)^ 
^- (l + 2£)(3 + 2^) ' ^^■^^> 

which yields G = 6/5, 48/35, and 100/63 for £ = 1, 2, 3. In order to see whether the 
cosmic variance or the intrinsic ellipticity distribution dominates the noise, we consider 
the ratio 

_ 2(Mfp)iV _ 15007r , f 9 ^ / a, x / n \ 
P - alG - -G- ^^^P^^^> 1,11^ ) Uoarcmin-^ ) ' ^^'^^^ 

where n = N/{n9'^) is the mean density of galaxy images. If p ^ 1, the intrinsic 
ellipticity distribution contributes mostly to the noise, whereas in the other case, the 
cosmic variance is the dominating factor. In this discussion the kurtosis term in (5.16) 
has been ignored. 

The coefficient p is plotted in the upper panel of Fig. 7, as a function of 9, for the 
same cosmological models as considered before. One sees that for filter scales exceeding a 
few arcminutes, the intrinsic ellipticity distribution of galaxies is no longer the dominant 
source of noise, but the 'cosmic variance' will start to dominate. 

5.2 Dispersion for an ensemble of fields 
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Fig. 7. The upper panel displays the quantity p - see (5.18) - as a function of angle 6, which describes 
the relative importance of the cosmic variance relative to the intrinsic galaxy ellipticity dispersion, for the 
same cosmological models as shown in Fig. 1. In this figure, we have set <Te = 0.2 and n = 30arcmin~^; 
the dependence on these two quantities can be seen from (5.18). Also, we have used the filter with I — 1. 
For all cosmological models considered here, the cosmic shear dominates for angular scales above a few 
arcminutes. In the lower panel, we plot the signal-to-noise ratio (5.22) for the dispersion on angular scale 
0, for the case that one has a square field of length L = 1 degree available, and that this field is densely 
covered with circles of radius 9, so that TVf = (L/20)^. We have assumed zero kurtosis for this plot. The 
signal-to-noise ratio at fixed is proportional to L. Parameters were chosen as in the upper panel. For 
both panels, the redshift distribution was chosen according to (2.18), with zq = 1 and j3 = 1.5, and the 
fully non-linear power spectrum was used 
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Obviously, from observing a single field with radius 6*, no reliable estimate for (M^p(^)) 
can be obtained. We next consider a sample of A^f fields, and assume that they are spa- 
tially sufficiently well separated so that the shear in one field is statistically independent 
of that in the others. In that case, an unbiased estimator of {M^^j is provided by 



A^= J]an 5^a,M, , (5.19) 



\n=l / n=l 

where is the estimator (5.2) in a single field, and the are weight factors. It is 
easy to show that the dispersion of M. is minimized if oc cr~^(M„), in which case it 
becomes 

\n=l / n=l 

If all A'f fields contain the same number of galaxies, then, as expected. 



^. (5.21, 



Combining (5.16) with (5.18) and (5.21), we find for the signal-to-noise ratio for a mea- 
surement of (Mfp): 

« - - ^ (5.22) 



N ■ a{M) 



//4 + 2(l + i)' 



The question of how far two fields have to be separated before they can be considered 
statistically independent can be investigated by considering the correlation between the 
values of Map in two fields of angular radius 9i and 6*2, separated by A9. We consider 
the correlation coefficient 

rM..Ae):=^-MsEMhmm, (5.23) 

'(M|p(9i)) {Mi^ie,}) 

which can be easily calculated with (3.5) and (3.2) to yield 

rcorri9i,92;A9) = . / ds s P,{s) Jo{s A9) h{s9^) h{s92) . 

■I{MU0,)){MI^{92))J^ 



(5.24) 

This can be compared with the correlation of the mean shear '^{9) as defined in (3.16), 
for which the corresponding correlation coefficient is the same as in (5.24) with Ii{sR) 
replaced by Itii{s9). In Fig. 8 we have plotted rcorr for various values of 6*1, 6*2 and A9, 
both for the Map-statistics and for 7. From the figure is can be easily seen that the Map- 
statistics decorrelates very quickly. For example, considering two circles of equal radius, 
such that they just touch (so that 9i — 92 — A9/2), we infer from the upper left panel 
in Fig. 8 that the correlation between the Map measurements in these two apertures is 
less than 1%! Therefore, if we had a large image, we can place apertures densely on that 
image and consider the Map-values obtained from each aperture as independent. Also, 
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for different values of 9i and 6*2, the decorrelation is very quick. This property can then 
be employed for obtaining (^M^^{9)'^ from a big image for various angular scales. 

Thus, for an image with sidelength L, we can place A^f = {L/29)'^ nearly independent 
apertures on this field. Using this estimate, we have plotted S/N, as defined in (5.22), 
in the lower panel of Fig. 7, assuming L = 1 degree. There we can see that the S/N ratio 
is larger than one for angular scales below 20'. For different values of L, one uses the 
scaling S/N oc L. The results of Fig. 7 can be compared with the Fig. 10 of JS which 
gives the signal-to-noise of the rms shear with a top-hat window. 

In contrast to Map, 7 decorrelates very slowly, as indicated by the light curves in 
Fig. 8. This can be traced directly to the shape of the top-hat filter Ith, displayed in 
Fig. 2, which picks up long-scale power of the density field. Therefore, if one wants to 
measure cosmic shear using 7, one has to use data fields which are well separated on the 
sky. 

5.3 Alternative estimators 

An apparently unrelated estimator for (M^p) can be obtained from the following consid- 
eration: define 

7(^) := J d^^ C/ (1^)7(1?) , (5.25) 

where the integral extend over a circle with radius 9. This definition is analogous to the 
one in (3.16), except that a weight function is added. Then 

(77*) = J d^9' U{9') j d^d U{d) (7(6>')7*(t?)> . (5.26) 

Noting that the two-point correlation function of the shear is the same as the two-point 
correlation function of n (see, e.g., Blandford et al. 1991), we see that this expression 
agrees with (3.10), so that 

{muo)) = mne)) ■ (5.27) 



Hence, a practical estimator for (M^p) is 

^ ^ NiN - 1) ^ ^^^'^ ^^^'^ ^''^''^ + ''^'''^ ■ ^^-^^^ 



With calculation similar to those in Sect. 5.1, one can show that indeed M is an unbiased 
estimator for (^M^^y, its dispersion in the same approximation as in (5.16) is 

a^(M)=A4(Mfp>V(^^ + V2(Mfp>j , (5.29) 

where p,4 is the kurtosis of the shear 7, and 

G = 7r9^ J d^^ U'^i'd) = G , (5.30) 

where the final equality is valid if U is chosen as in (3.12). Thus we see that the dispersion 
of M is very similar to that of M, except for a factor in the term containing a^- 
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Fig. 8. The correlation coefficient rcorr(^i, ^2; ^d) is plotted as a function of A9/9i (heavy curves). Here, 
61 is the larger of the two aperture sizes; the other one is 1, 1/2, 1/4, and 1/8 of the former one in the 
four different panels. Six curves are drawn in each panel, corresponding to 0i/arcmin = 1 (solid curves), 
2 (dotted curves), 4 (short-dashed curves), 8 (long-dashed curves), 16 (short-dashed-dotted curves), and 
32 (long-dashed-dotted curves) . The light curves in each panel is the corresponding correlation coefficient 
for the mean shear inside circles, y{9). For this figure, an i?^ — 1, i7v — universe has been assumed, 
with as = 0.6 and F = 0.25, and the sources were assumed to follow the redshift distribution (2.18) with 
zq = 1 and P = 1.5 



This difference is due to tfie fact tliat in M, two components of e are used, each of 
which carries its dispersion. Hence, for measuring (M^p(6')), one can use the estimators 
(5.2) and (5.28), where the former should yield a slightly higher signal-to-noise ratio. A 
comparison between these two estimates can be used to check the integrity of the data 
and the data analysis procedure. 

Another check on the quality of the data can be obtained by noting that 
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Mr:= J d2^Q(|i?|)7,(t?) (5.31) 



should vanish identically, where 7^ is the radial component of the shear, defined by taking 
the imaginary part in (3.9) instead of the real part for the tangential component. The 
fact that Mr = can be shown easily by introducing the Fourier transform of 7r to obtain 
a result similar to (A5), with the cosine replaced by a sine; the polar angle integral in 
(5.31) then yields zero. The estimator 




should therefore yield a value within the dispersion, which is given by a{M.j-) = 
alG/{^/2N'iN). 

5.4 Practical estimator for (M^p) 

An obvious estimator for (M^p(6')) for a single field is 

Ms := QiQjQkeuetjetk ■ (5.32) 

i,j,k=l 

It is easy to show, using the methods used before, that E{Ms) = (Mfp). An estimator 
for an ensemble of fields can then be defined immediately in analogy to (5.19). 



6 Discussion and conclusions 

In this paper we have introduced the Map-statistics, or mass-aperture statistics, as a 
new measure for cosmic shear. We have compared those statistics with others suggested 
earlier, namely the shear two-point correlation function, and the rms of the shear averaged 
over circles. The main results can be summarized as follows: 

(1) Map is defined as a filtered version of the projected density field k, but can be 
calculated in terms of the tangential shear inside a circular aperture. Hence, Map has a 
well-defined physical interpretation, and at the same time can be expressed in terms of 
the observable shear. 

(2) The dispersion of Map, (Mfp(6')), can be expressed as a convolution of the power 
spectrum of the projected density field with a filter function which is strongly peaked. In 
contrast, the corresponding filter function for the rms of the shear averaged over circles is 
considerably broader and in particular picks up the long wavelength range of the power 
spectrum. Therefore, the Map-statistics has the nice property of being a well localized 
measure for the projected density field. From measurements on different angular scales, 
the power spectrum of the projected density field can be constructed. A different method 
to obtain a local estimate of the power spectrum of the projected density field has recently 
been suggested by Kaiser (1996). 

(3) In contrast to the mean shear inside a circle. Map is a scalar quantity. It is 
therefore possible to define the skewness of Map which is measurable. 
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(4) Whereas (Mfp(6')) is smaller than ^|7(^)|^y, at least on small scales, the dis- 
persion of the estimator for (M^p) is smaller by a factor of \f2 than that for ^|7(^)|^^ 

on small scales, due to the fact that in the former case, only one component of the shear 
is needed for the estimate. On scales beyond a few arcminutes, the dispersion of the 
estimates for both statistics is dominated by the cosmic variance, so that the fractional 

— 1/2 

accuracy of the estimates is nearly ' , where A^f is the number of independent circular 
apertures. 

(5) The values of Map for neighboring circles dccorrclates very rapidly with increasing 
separation. In particular, we found that the values of Map calculated on circles which just 
touch are nearly mutually independent, with a correlation coefficient below 10~^. This 
implies that one can make use of wide-field images by 'punching' circles on these images 
and considering the Map values on these circles as mutually independent. In contrast to 
that, the mean shear inside neighboring circles is very strongly correlated, and these have 
to be separated by several degrees before they can be considered independent. Therefore, 
the Map-statistics provides a very valuable tool for measuring cosmic shear on scales of 
arcminutes from wide-field images. A single one square degree image of sufficient depth 
(e.g., corresponding to a few hours integration time on a 4-meter-class telescope) and 
high imaging quality can be used to obtain (M^p(6')) on scales below ~ 5' with a signal- 
to-noise larger than ~ 3. The signal-to-noise is proportional to the square root of the 
solid angle covered by the available data. 

(6) The skewness of Map is very sensitive to the cosmological model, but rather 
independent of the shape of the power spectrum. In particular, in the quasilinear regime, 
the skewness is independent of the normalization of the power spectrum. Whereas the 
latter property may not be strictly preserved when considering the skewness in a fully 
non-linearly evolved density field, the experience from comparing numerical simulations 
of the density field with quasi-linear predictions has shown that at least the skewness 
of the density field is rather well approximated by quasi-linear theory, even for density 
contrasts of order unity (Baugh et al 1995, Colombi et al. 1996, Gaztaiiaga & Bernardeau 
1997). We therefore expect that the independence of the skewness on the normalization 
will be roughly preserved in the non-linear case. This makes the skewness a very valuable 
probe for cosmological parameters. 

(7) The difference between the predictions from linear evolution of the cosmic density 

fluctuations, and the fully non-linear evolution of the power spectrum, is much larger for 
the m-statistics than for the mean shear within circles - see Fig. 3. This is related to 
the fact that the filter function for (Mfp) is narrowly peaked, whereas the corresponding 

filter for ^|7(^)|^^ has a long tail towards long wavelengths. Since the long wavelengths 

have not become non-linear (see Fig. 1), the non-linear effects for ^|7(^)|^^ are necessarily 

weaker than for (M^^p). 

The signal-to-noise estimates given in this paper suffer from the lack of knowledge 
on the kurtosis of Map, which we have not attempted to calculate, or more generally, the 
distribution function of Map. Therefore, the estimate presented in Fig. 7 may be slightly 
optimistic. A more accurate estimate can be obtained from numerical simulations, by 
taking high-resolution realizations of the cosmic mass distribution and studying light 
propagation through such a Universe. Such studies have been undertaken in various ways 
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and with various scientific goals in the past (e.g., Jaroszyhski et al. 1990, Lee & Paczyhski 
1990, Bartelmann & Schneider 1991, Jaroszyhski 1991, 1992, Wambsganss et al. 1995, 
1997). In particular, Blandford et al. (1991) have compared their analytical estimates 
for the rms shear inside circles with results from numerical ray propagation simulations. 
Further work in this direction will most useful. In particular, the performance of practical 
estimators for the dispersions and skewness can be evaluated, and the accuracy of the 
approximations used here (i.e., the Born approximation and the neglect of lens-lens terms 
in (2.10)) can be checked. 

In summary, the Map-statistics appears to be an attractive measure for cosmic shear. 
We expect that its applications can be expanded beyond the range considered here. It 
is a promising method to search for mass concentrations in the Universe, as pointed 
out in Schneider (1996). A further application may include the correlation of shear, as 
measured by Map, with the distribution of foreground galaxies, measured through the 
same filter function U. We shall exploit this route as a method for measuring the bias 
parameter, and its dependence on scale, in a later publication. 
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for a careful reading of the manuscript. This work was supported by the "Sonder- 
forschungsbereich 375-95 fiir Astro-Teilchenphysik" der Deutschen Forschungsgemein- 
schaft, by the Programme National de Cosmologie of the Centre National de la Recherche 
Scientifique in France, and by EU contract CHRX CT930120. 



Appendix 

In this Appendix we calculate the contributions to the (observable) skewness which would 
be present even in the case of a Gaussian density field. The three effects which cause 
a non-zero value for (M^p) also for Gaussian fields are: (1) Image ellipticities provide 
an unbiased estimate of the reduced shear g — 7/(1 — re), rather than the shear 7 
itself. (2) The separation vector x(0,w) deviates from fK{w)6 - see (2.7); the ('Born') 
approximation leading from (2.9) to (2.10) neglects this effect. (3) Also, in the same 
step, the matrix A on the right-hand-side of (2.9) was approximated by the unit matrix; 
therefore, (2.10) does not contain the coupling between defiectors at different redshifts. 
The two latter terms have already been discussed, in somewhat different context, in 
BvWM. We shall consider here these effects in turn. In this Appendix, 5 and are 
meant to be the linear density field and its corresponding gravitational potential, since 
we are interested in the contribution to the skewness coming from the above-mentioned 
effects in the presence of a Gaussian field. 



Al Aperture mass in terms of g 

Let 

be the tangential component of the reduced shear. Then, from image ellipticities, one 
obtains an unbiased estimate of 
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where 

5M,^{9) ^ j g(|t?|) 7t W • (^3) 

Hence, the observable skewness becomes 

(M|(^)> « (M3p(^)> + 3 (Mfp(^) <JMap(^)> ^ {M%{e)) + <5 (M,^p(^)>^ . (yl4) 

Using the fact that the Fourier transform of the shear is 7(3) = k{s) e^^*^", where is 
the polar angle of s, the tangential shear - see (3.9) - can be Fourier-decomposed as 

7tW = -7^e(| ^e^^-=«(s)e'«-^e-^^^) =-|^cos2(^-^,)e^«-^«(s), 

(^5) 

where in the last step use was made of the fact that ^(1?) is real, and (f is the polar 
angle of "d. Then, from the definition in (A4) and by introducing the expression (3.5) for 
Map (6*), one obtains after replacing k('i?) by its Fourier representation: 



X 



I d^^Q(|t?|) j |^e'«'-^y cos2(<^-<^,)e'«-^(«(sO^(s2)«(s)K(s')) • 

(A6) 

We calculate the correction 5 (Mfp)^ to leading order and thus use the linear evolution of 

the cosmic density fluctuations. In this approximation the density field remains Gaussian, 
and since k is a linear functional of the density fluctuations, k is Gaussian. Then the 
correlator in (A6) becomes 

(«(si)«(s2)«(s)«(s')) = (27r)^P«(si) P«(s') 5d(si + S2) 5d(s' + s) ^^^^ 

+ 2 terms obtained from permutation . 

Inserting (A7) into (A6), one notices that the first term yields zero, whereas the other 
two terms yield equal contributions. After carrying out the Si and S2 integrations, one 
finds 



X j cos2(v.-v..)e^-(^-^^)P.(.) / (|;^e'^'-(^-^^)p.(.') . 



Next, we perform the polar angle integration of i?i and ■i?2, and use the definition of the 
filter I (,{7]) to obtain 



5(MiW>, = -6 / d^^g(|t?|) j -^cos2{^-^s)^^-'^P.{s)h{se) 
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Finally, the three integrations over the polar angles can be carried out, leaving the final 
result 



HM!p{0)}g = 127r j dssP^{s)Ie{se) j ds' s' P,{s') h{s'e) 
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Fig. 9. The fractional change of the skewness which occurs by considering ^Mfp^ instead of {|m|^, 
where the latter one is the observable quantity, for four different cosmological models (thick curves). 
For comparison, the thin curves show the corrections due to the Born approximation and the neglect 
of lens-lens coupling terms in (2.10). In this figure, the sources were assumed to follow the redshift 
distribution (2.18), with 2:0 = 1 and (3 = 1.5 




In Fig. 9 we have plotted the fractional contribution (Mfp(6'))^ / (M^p(6')) as a 
function of angular scale 6, for four different cosmological models. As can be seen, the 
difference between (M^) and (Af^p) is very small, lower than 6% for all cases considered, 
where the largest deviations occur in the EdS models. 

A2 Lens-lens coupling, and dropping the Born approximation 

The expression (2.10) for the Jacobi matrix A is valid up to first order in the gravitational 
potential The two effects considered here are obtained by expanding (2.9) up to second 
order in 

We write (2.7) as x(6>, w) = x(o) (6>, «;)+x(i) (6>, w)+0{$'^), with x^") (6>, w) = /k(w)6> 
andx(^)(6>,w;) given by the second term in (2.7), now with x in the integrand replaced by 
x(0). Similarily, we write A{e, w) = ^(°)(6>, w) + A^^\e, w) + A^^\e, w) + 0{^^), with 
A\f{0,w)=d,j, and A^^\0,w) given by the second term of (2.10). Expanding (2.9) in 
terms of <5, one finds that 
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dw 



f fK{w~ w')fK{w') 



X 



/^H (^11) 

^^iki {fK{w')d,w') x{^\d,w')5kj + ^,ik {fK{w')d,w')A]^j^{d,w')] , 



which, after inserting the expressions for x^-*^) and A^^^ yields 



4 

c Jo 



, fxiw- w')fK{w') 



X 



+ 



fK{w) 

UK{w')e,w') ^,1 {fK{w")e,w") 

fK{w") 



i-w' 

/ dw" fK{w' — w") 
Jo 



iAl2) 



i^,ik {fK{w')d,w')^,kj {fK{w")e,w") 



We consider 1/2 of the trace of A^"^^ as the surface mass density to second order in 
Although A^'^^ is in general not symmetric, so that to this order the lens mapping can no 

longer be described by an equivalent surface mass density, the asymmetry is expected to 
be small in realistic situations, and we shall neglect it henceforth. Then, after averaging 
over a source redshift distribution, as done in (2.15), and replacing the derivatives of $ 
by their Fourier representation, using the Poisson equation (2.13), one finds 



dw 



„ fK{w'-w") 



a{w") 



X 



X 



X 



d^fc' 1 



exp [i {kW + fK{w')\^ ■ e)] 5{k'- w') 

d^k" 1 



(27r)3 |^//|2 



exp [i {k'iw" + fK{w")\^' ■ d)] 6{k"; w") 



M!f!;i(k'.k'f + |k'|^k'.k' 

fK{w') 



{A13) 



where we use the same notation as in the main text, i.e., k has k as the first two 
components. 

To calculate the resulting contribution to the skewness according to (4.1), we note 

that 



(AU) 



where the '='-sign means that the two expressions yield the same result after insertion 
into (4.1), owing to the symmetry of the integration there. The correlator in (A14) is 
calculated by inserting (A13) for k^'^\ and (2.15), with 5 replaced by its Fourier decom- 
position, for K^^\ The resulting expression then contains a correlation function of four 
5-terms. Since we consider the lowest, i.e., the linear order of S, the four-point correlation 
function is given by (4.9). Using arguments very similar to those which led from (4.3) to 
(4.5), one finds 
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a^(t(;') 



X 



^ Po{k') exp (i/K(«;')k' • (02 - 03)) 



Mw) Ik^k'!' |kf , 



(A15) 



As a final step, we change variables to s = /K(w)k, s' = fK{w')k.' , and insert the resulting 
expression into (4.1), using (A14). After applying the definition of Ii, the contribution 
to the skewness from lens-lens coupling and dropping the Born approximation becomes 



r^^, g{w')fK{w~w')D\{w') 
Jo a^{w')fK{w) 



X j d^s Po 



X 



jKiw) 



(s ■ s')^ s ■ s' 

I |2 I ,|2 I ,|2 



(^16) 



./k(«'0 



We have plotted this correction term, together with 5(Mfp)^, in Fig. 9. There we can 
see that this correction term is of the same order as that considered in the previous 
subsection, i.e., smaller than ~ 5% on scales larger than 1 arcmin. These corrections 
are therefore smaller than the uncertainties introduced by calculating the skewness with 
quasi-linear theory. A numerical ray-trace calculation would of course take all the cor- 
rection effects mentioned here into account. 

It turns out that AS is nearly independent of cosmology; this can be traced back to 
the fact that d (Mfp) and (Mfp) can be expressed as bilinear functions of the projected 
power spectrum P^. Assume for a moment that locally, Pk{s) is a power law; then AS 
would depend only on the local slope of this power law, independent of cosmological 
factors. Hence, if the local slopes of the power spectra in different cosmological models 
are similar, one expects AS to be nearly independent of the cosmological model. 
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